LOW-MACH-NUMBER EULER EQUATIONS WITH SOLID- WALL 
BOUNDARY CONDITION AND GENERAL INITIAL DATA 



BIN CHENG 

Abstract. We prove that the divergence-free component of the compressible Euler equations 
with solid-wall boundary condition converges strongly towards the incompressible Euler equa- 
tions at the same order as the Mach number. General initial data are considered and are not 
necessarily close to the divergence-free state. Thus, large amplitude of fast oscillations per- 
sist and interact through nonlinear coupling without any dissipative or dispersive machanism. 
It is then shown, however, that the contribution from fast oscillations to the slow dynamics 
through nonlinear coupling is of the same order as the Mach number when averaged in time. 
The structural condition of a vorticity equation plays a key role in our argument. 



1. Introduction 

Hyperbolic Partial Differential Equations (PDEs) of multiscale nature have seen rapidly grow- 
ing applications in recent years. We conduct a theoretical investigation in this paper concerning 
a prototypical example: the compressible Euler equations of barotropic fluids in a bounded do- 
main in the low-Mach-number regime. The system of equations, in terms of total density p and 
velocity u, read 

d t p + V-(pu) = 0, 

p(d t u + u-Vu) + ^m^ = Q, 
where Mach number e<l brings in fast scale of oscillations. The equation of state is 
(1.1) pressure=p(p) with p(-) £ C°°, p(l) = 0, p'(l) = 1. 

The spatial domain Q, £ R D (D = 2 or 3) is bounded with smooth boundary <9f2. It is connected 
but not necessarily simply connected. The solid-wall boundary condition is prescribed as 

Hen = °- 

Here and below, n = n(x) denotes the outward normal at x £ dil. Without loss of generality, 
impose 

1 f , 
— — / pax = 1 



n 



\n 

since f^pdx is always conserved by the dynamics. 
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Having p ~ 1 as the total density, define the density perturbation 

P — 1 

(1.2) p:= P — , 

and rewrite the above system for the unknown pair (p, u), 



(1.3) 



V-u 

d t p + V-(pu) + — = 

„ p'(l + ep) Vp „ 

d t u + tt-Vtt + — — — = 

1 + ep e 



(1.4) u-n\ dQ = , / pdx = 0. 



Theorem 1.1. Consider the initial-boundary value problem of the D- dimensional compressible 
Euler equations (jl.3p . (jl.4p subject to initial data (po, uq) 6 H m {yt) with m > y + 4. Assume 
(Po,uq) is compatible with the boundary condition "d^uo"-n\g n = /or k < m. Here, "d^uo" 
is obtained from expressing d^u in terms of spatial derivatives via (|1.3p and i/ien substituting 
(po, no) into i/te expression. 

Then, there exist general constants E, T, C only dependent on m, $7 and £/ie pressure law 
p() s.t. for all e £ (0, 2£/||(po, ^o)||ff m ]> ^ fte solution (p,u) exists smoothly for times t £ 
[0,r/||(po, no)||j3"™] and ||(p, u)||ffm < C||(po, ^o)||i? m uniformly for such times. 

Moreover, there exists a splitting of the velocity field u = u p + with V-u p = and = Vtjj 
for some potential function vp s.t. for general constants c\,C2 only dependent on m,Q,p(-) and 
||(po, uo)\\h™, 

• (strong convergence of the slow compoenent ) 

o maxJn P (t,-) - u(t, -)||ff™-3(n) < c\E 



(1.5) 



where u solves the incompressible Euler equations 

{ d t u + u-Vu + Vg = 0, 



u 

and, 



V-u = 0, 
p 



n laa = ' u(0,-) = u 



• (weak* convergence of the acoustic part) for any smooth testing function (p' , u') defined 
on [0, T] x Ti, 

(1.6) ft pp' + u Q -u'dxdt\ < c 2 e(\\(p', u')\\ LTLl + \\d t (p', u')\\ L 2 ). 

io in I 

Proof. For strong convergence, see Theorem 13.11 Some major new ideas of this article are 
centered around this theorem. For weak* convergence, see Corollary 12.61 

For uniform existence time and H m estimates, see Theorem 14.61 and the remarks thereafter. 
Note that the mere existence of H m solution is not the concern of this article. Thus, the 
compatibility condition serves an assumption in this theorem and will not appear explicitly 
again. □ 
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There have been numerous results regarding the singular limits of compressible Euler equa- 
tions and other fluid equations in various settings. We point to two survey papers for some 
comprehensive lists of references: Schochet [23] with emphases on hyperbolic PDEs and homog- 
enization in space-time; Masmoudi [17] with emphases on viscous fluids and weak solutions. To 
mention only a few earliest work, we note papers by Ebin [TJ [9] , Kreiss [14] , Beirao de Veiga [2] , 
and Klainerman and Majda [13]. In a closely related paper [JJ, Kreiss etc. applied the bounded 
derivative method in numerical schemes to gain control on time derivatives and thus to rid of 
fast gravity waves. These results, in terms of (|1.3j) . confirmed that compressible flows (p, u) 
converges to (0, u) strongly at order 0(e) with u solving (jl.5p provided the initial data (po, uq) 
also converges to (0, uq) strongly at order 0(e) 

(1.7) \\po\\ + 1 1 i*o — || ~ £ for some div-free uq. 

Here, || • || denotes suitable spatial norm. Note that condition (|1.7p implies that perturbation in 
the total density p vanishes at order 0(e 2 ) — consult (jl.2j) . 

This family of well-prepared initial data lead to uniform estimates on the size of dt(p, u) at 
t = 0, independent of e, by the virtue of fjl .3|) and therefore suppresses the so-called initial layer. 
Then, one obtains uniform control on the size of dt(p, u) for finite times, which allows passing 
of limits by the Arzela-Ascoli Theorem. Well-prepared conditions on initial data were later 
removed for problems in the whole space (Ukai |25j), in an exterior domain (Isozaki [10\ 111] ) 
and in a torus (Schochet [22] ). These arguments more or less rely on use of Fourier analysis 
and/or dispersive nature of the underlying wave equations. 

Singular limit problems in a bounded spatial domain, on the other hand, remain much less 
studied. Schochet proved in [20J the same low-Mach-number limit with solid-wall boundary 
condition and, again, well-prepared initial data. A main challenge in this setting is the presence 
of characteristic boundary. It is elaborated in Rauch's work [18] for linear systems that, in 
general, only estimates along tangential directions are available near the boundary. We also 
note that there were also preceding results in e.g. [9], [2], all of which required well-prepared 
initial data. 

The main originality of our paper is to prove the low-Mach-number limit for general initial 
data without requiring the well-prepared condition (II. 7h . We establish a strong convergence rate 
of 0(e) for the slow incompressible part and weak* convergence for the fast acoustic part. In 
lieu of the commonly accepted initial layer, it is proved that nonlinear resonance of fast acoustic 
waves does not enter the slow dynamics at all (Lemma l3.3p and interaction of fast-slow dynamics 
vanish at order 0(e) upon integrating in time (Lemma I3.4p . This argument was explicitly used 
in Cheng [5] to study the rapidly rotating shallow water equations with 2 fast scales in the whole 
space. It recently came to our knowledge that similar ideas have occasionally appeared in the 
literature for weak limits (cf. equations after (7) in Lions and Masmoudi [16]) and for problems 
without boundary (cf. equation (4.27), (4.28) in Schochet [23]). We also tackle the boundary 
condition carefully and present a clear calculation of a priori estimates on space-time norms of 
the solution. Here, the vorticity equation plays a crucial role, which was argued in e.g. Schochet 
|21j . Throughout our analysis (not just a priori estimates), we employ elliptic estimates for 
PDE systems with general boundary conditions (Agmon etc. [1]). 
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The organization of the rest of this article is as follows. In Section [21 we introduce elliptic esti- 
mates particularly for the case at hand and present a precise characterization for the projection 
operators associated with Helmholtz decomposition. This projection is to be used to split the 
solution as well as the system into "slow" and "fast" parts. Moving on to Section [3j we conduct a 
thorough study on nonlinear interactions of fast-fast and fast-slow types, which will confirm the 
assertions in the previous paragraph. Solution regularity is presumed here and is then treated 
rigorously in the next Section H] regarding a priori estimates. A mollification method is used 
to handle the lack of boundary regularity. The proofs are self-contained, only relying on basic 
Calculus. Finally, Section [5] contains a crash course on the rapidly Rotating Shallow Water 
equations to which the exact same ideas apply. It is presented in a very systematic way due to 
an algebraic structure of duality nature. 

We will repeatedly use some well-known inequalities of Sobolev norms without making refer- 
ences. They are all based on Holder's inequality, Gargliardo-Nirenberg inequality and Sobolev 
inequality. For the most part, it is sufficient to accept the following estimates, 

(1-8) \mgdl 2 g...di*g\\ L2{Q) < c\\g\\ 

where m > D/2 + 1, < j\ < ... < jk < m and ji + ... + jk < m + 1. 



2. Elliptic Estimates and Helmholtz Decomposition 

Elliptic estimates and Helmholtz decomposition are both from Elliptic PDE theories, the 
former regarding regularity and the latter solvability. In connection with the compressible Euler 
equations, the singular terms in (|1 .3j) with fast scale \je define an elliptic operator, 



C 



p\ = fV-u 
u) ' \Vp 



and there is another elliptic operator 



I,] I P \ ■= ( _° ] in :;D and I I in 2D 




u) \Vxu 

that essentially yields the vorticity. They satisfy, 

KC = and, in 3D, OC = 0, 

which is why the singular 1/e term does not appear in the vorticity equation. 

The papers of Agmon, Doughs and Nirenberg [lj establish a Complementing Boundary Con- 
dition that is necessary and sufficient for the solution operator of a s-th order elliptic PDE 
system to be C m C m+S and H m — > H m+S . In this article, only a particular case is used: for 
any velocity field u with a trace subject to the solid- wall boundary condition u-n|g^ = 0, 

(2- 1 ) IMk™(n) < C (HV-wllij-m-i^) + ||Vxu||#m-i( n) + ||n|| L 2 {Q) ) . 

Here and below, we always assume m is a positive integer so that the trace u\ d ^ is well-defined. 
See e.g. [3] for application of this estimate. 
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Remark 2.1. The ||it||/,2 term in the above estimate is removable if and only if the elliptic PDE 

system 

(2.2) Vxu = 0, V-u = u-n\ m = 0, 

admit a unique trivial solution u = 0. This is the case e.g. if f2 is contractible (in virtue of 
the Poincare lemma) but is not true in general. For example, in 2D, when d£l consists at least 
two disjoint branches Ti, T2, one solves a scalar Laplace equation = with ip\ =04 and 

a\ ^ ci2- Then, u = ( J is indeed a nontrivial solution to the above PDE syste We will 
\-d x ipJ 

revisit this remark two more times below. 

For the p component of the solution, under the zero mean condition in (jl.4p . one has 
Hp|lL 2 (f2) ^ CII^pIIl 2 ^) by the Poincare inequality and therefore 

( 2 -3) ||p||h™(Q) < C\\Vp\\ Hm -i( Q y 

Now, define a solution space X m C H m (Q) as 

X m := UJ = Q G H m (Q) u-n = on dtl and J pdx = j . 



Then, the above estimates (|2.ip . (|2.3|) lead to 

(2.4) \\U\\ Hm{fl) < C(\\CU\\ H m- 1{n) + II^H^-i^j + ||n|| L2(Q) ) for U € X™ 

Next, we move on to Helmholtz decomposition. Define V as the 1? projection onto the I? 
closure of Ker£ n X m and define Q as its orthogonal complement, 

(2.5) V := L 2 -Proj{Ker£nX mL2 }, Q-=Z-V. 

Note that Ker£ n X m = | \®\ G H m u-n\ QQ = 0, V-u = 

We present two ways to characterize V, Q. The first way is to express them using an elliptic 
PDE, 

(2.6) for any C7 = (f) 6 JT-, QV = (> 



with <f> solving 
(2.7) 



Aft = V-u, 
V4>-n\ m = u-n\ dn . 



Proposition 2.2. Operator Q given in (|2.6p . (|2.7p exists uniquely and amounts toZ — V when 
restricted on H m (£l). 



With a compactness argument, one can further show that these solutions form a finite dimension space. Then, 
(|2.1[) can be expressed in quotient space with the L 2 norm removed from the RHS. 
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Proof. It is sufficient to study the well-posedness of (|2.7p . Define a smooth function 

4>i(x) := x(dist(x, dQ))u(x) -Vdist(x, 8ft) for 

where %(•) G C°° s.t. with a suitable constant c, x( z ) = 1 if |#| < c and x( z ) = if |z[ > 2c. 
It is easy to verify that V^i-n]^ = li-n]^ since Vdist(x, = — n. Then, one uses the 

Fredholm alternative to show that 

( A(j) 2 = V-u - A</>i 
\V0 2 -n|^ = O 

admits a H m+l solution 4>2 that is unique module a constant. In fact, this is done by showing 

/ (Vu - Afa) dx = (u-n - V0i-n) (is = 0. Therefore, V</> = V(0i + 2 ) in (|2"1)J) . (I2T71) exists 
Jon 

uniquely in the space of 

Next, we verify Q = I — V. By definition, 

and V(u — V(/>)-n| 9n = 0. 

Therefore, 

(I - Q)C/ e Ker£ n X m . 
It remains to show QU is orthogonal to ImgV. Take a testing function U' = \ ,) € Ker£nX r 



that satisfies p' = and V-i/ = with ^ /- n|gQ = 0. Combining this information, we calculate 
/ U'-QUdx = / p p + u'-V(j)dx = \ V-(u'4>)dx= / wu'<j)ds = 0. 

Jn Jn 
Thus, Qf/ is orthogonal to Ker£ n X m and therefore its I? closure. □ 

Now that V, Q are well-defined, for simplicity, we will use U p for VU and U® for QU 
whenever it is not ambiguous. We will also use lower case u p , ifi for the associated velocity 
components. The density component of VU is always zero. 

The second way of characterizing V essentially relies on the skew-self-adjointness of C and 
will be very useful in studying nonlinear resonance of acoustic waves. 

Proposition 2.3. In addition to the definition ofV, a duality relation also holds true, 

(2.8) KerP = Img£ 

JJm X m + 1 

Proof. Since V is a projection, its kernel identifies with Img(Z — V). So it suffices to show 



ImgQ = Img£ 



Take any U Q = (^\ G H m . By (22)]), ([ZZ3), u Q = V<f> for some G H m+1 with J Q $ = 0. 
Also, since f n p = 0, the Poisson's equation Aip = p with V^-n|^ = admits (at least) a 
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solution in H m+2 . Thus, we set U' = ( „ ,) G X m+1 and clearly U Q = CU'. Therefore, 



LHS(ES)CRHS(ES 



Assume U = J = CU' for some U' = K,J G X m+1 . Since u = Vp', by (J2SI), ([23D and its 

unique solvability, we have u Q = it. Thus, U = U Q G ImgQ. Therefore, RHS(|2^cLHS([lLH]). 

□ 

Remark 2.4. The above duality relation reminds us of Linear Algebra if C, K, were (skew)- 
symmetric matrices and V, Q were matrices representing orthogonal projection onto their null 
spaces. In fact, if is contractible, one can use the Poincare Lemma to show that KerV = 
Img£ = Ker/C in proper regularity spaces. In general, however, only KerP = Img£ C Ker/C is 
true. Take for example u given in (|'2.'2p . It is a nonzero element of both Ker/C and KeiC = ImgP 
and therefore not of KerV. 

On the other hand, regardless of domain topology, both Img/3 and Ker/C are invariant spaces of 
the barotropic compressible Euler equations. In other words, if the initial velocity is a potential 
flow (resp. curl-free), then it stays so at later times for as long as the classical solution exists. 

Now, we claim that V, Q are both bounded operators in H m (Q). Indeed, by definition CV = 
and thus CQU = CU. Also, by Proposition 12.31 ImgQ = Img£ C Ker/C and thus JCQU = 0. 
Then, elliptic estimate (I2.4h implies 

(2.9) \\QU\\ H m < C(||£E/|| H m-l + \\QU\\ L 2) < C\\U\\ H m 

(2.10) and likewise \\VU\\ H ™ < C||f7||fl*». 
Even more can be said about Q. 

Proposition 2.5. For any U® G ImgQ n X' m (m > 1), 

\\U Q \\ H m < C\\CU Q \\ Hm -l. 

Proof. Inequality (|2.9p being established, it remains to estimate ||QC/"||z,2. 

By Proposition 12. 3( vP = (with <fi chosen to have zero mean). With the assumption 
Ian = 0' ^ follows 

f \u Q \ 2 = [ u Q -V<j)= [ V-u Q (j) 
Jn Jn Jn 

<\\^U Q \\ L 2U\\ L 2 <C||V-^|| L2 ||V^|| L 2. 

Since v9 = V(j), we deduce Hit^H^ < C|| V-u^H^a which obviously is bounded by C\\CU^\\ H m-i. 

□ 

A useful application of the above propositions is the weak* convergence result of Theorem 

o 



Corollary 2.6. Let C7 = ( ) solves the compressible Euler equations fll .3f) . (II. 4h with sufficient 

W 

regularity. Then, there exists constant C depending on the size of U s. t. for any testing function 



s 
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(2.11) 



/ / U Q -U'dxdt <Ce(\\U'\\ L oo L i + \\d t U'\\ L 2 ). 
Jo Jn ' x *' x 



Proof. Since we always have U^-VU' = 0, it is enough to consider U' = QU' . Then, by 
Proposition E31 there exists U" G C\ x n X 2 s.t. U' = CU". Since CU" = CQU" , choose 
U" G ImgQ as well. By Proposition 12. 5( 

(2-12) HfTlfli < C\\U'\\ L 2, \\d t U"\\ H i < C\\d t U'\\ L 2. 

Now, by the skew-self- adjointness of C w.r.t. to the solid-wall boundary condition, 

f U Q -U'dxdt= [ [ £U Q -U"dxdt= [ [ CU-U" dxdt. 
o Jn Jo Jn Jo Jn 



Then, use the compressible Euler equations (II. 3D to continue 

(d t U + bounded terms) -U" dxdt 



o Jn 



e / U-U"dx 



+ e 

o 



f [ (U-d t U" + bounded terms-C/") dxdt. 
Jo Jn 



Plug in (|2.12|) and the rest of the estimate follows. □ 

This weak* convergence result can be interpretated in terms of observations in practice. If 
one chooses V in (12. Ill) as an integral kernel 

U' = 4>(t — to, x — xq) with / (j)dxdt = l, 



and uses it to represent some kind of physical measurement of the velocity and density fields 
in the fluid, then (|2.1ip confirms that the presence of the fast, oscillatory component in the 
resulting measurement is of order 0(e). This is one of the reasons why (almost) incompressible 
fluid motions are often observed in practice. 

3. Estimates on Nonlinear Interaction and Strong Convergence 

Rewrite (11.31) in terms of U = ( 

(3.1) d t U + N{U,VU-e) = --C[U}, u-n| 9Q = 0, [p = ° 

£ Jn 

where the nonlinear term 

(? n\ A \f(TT V77T \ f MrVj02 + PlV-«2 

(3.2) and M{U^ W 2 ; e) := ^ + ^ ^ 

(3.3) with 9{p -e),= (^±^-^ 1 



1 + ep J e 

In this section, we always assume the solution to the above system exists and satisfies U G X r ' 
for some m > D/2 + 4. It then follows that d t U G X" 1 " 1 , / U dt G X m , etc. 
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Theorem 3.1. Consider the D- dimensional barotropic compressible Euler equations (|3.ip . As- 
sume a solution exists classically U(t,x) G n^gC ([0,T]; H m ~ 3 (fi)) with m > D/2 + 4. Also 
assume a non-vacuum condition ||ep||c([o,T];H m ) — 1/2- Then, there exists an incompressible 

flow U = (®\ s.t. for t G [0, T] 

\\PU - F|| m _ 3 < CeF(max ||C/|| m ), 

mt/i projection P defined in (12. 5p and some polynomial. 

In particular, u can be chosen as the unique H m solution to the incompressible Euler equations 



(3.4) 



dtu + u-Vu + Vq = 0, 
V-u = 0, 



u ' n \dn = > u(0,-) = Uq. 

In other words, U solves 



(3.5) dtU + VN(U, VU; 0) = 0. 

Remark 3.2. Here, we avoid using the vorticity formulation for the incompressible Euler equa- 
tions since vorticity and divergence alone do not necessarily determine a unique velocity field. 
Consult the example in (|2.2p . 

The key to proving this theorem is contained in two lemmas. They have been used in [S] in 
the torus domain for the rotating shallow water equations. We now migrate them to bounded 
domains using the properly defined projection operators V, Q. 

Apply V to d£TJ. % Proposition [231 VC = 0. Thus, 

d t U p + PAf(U, VC/; e) = 0. 

Although M defined in (j3.2j) is not entirely bilinear, we take advantage of the facts p'(l) = 1 
and p p = to derive 

d t U p +PN(U p , VU P ; 0) + PM(U Q , VU Q ;e) 

+PM{U Q ,VU p ;0) +PM{U p ,VU Q ;0) = 



In comparison with (13. 5p . there are two types of nonlinear interactions to be studied, 

"fast-fast" PM{U Q ,U Q ;e), 

"fast-slow" PM(U Q , U p ; 0) + PN{U P ', U Q ; 0). 

Intuitively, one expects the nonlinearity of M to generate resonance from all kinds of interac- 
tion among U p and U®. However, the following lemma excludes the contribution of "fast-fast" 
interaction from ImgP. 

Lemma 3.3. For any G ImgQ with sufficient regularity, 

PM{U Q ,VU Q ;e) = 0. 
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Proof. The density component of lmgV is always zero. By Proposition 12.31 h suffices to show 
that 

ift = V(j) implies vft-VvP + g(p; e)Vp is also a gradient. 
This is true due to the Calculus identity 

1. 



(V)0-V(V0) = -V|V0| 2 , 



and the fact g(p; e) is a function of p only. 



□ 



Next, move on to "fast-slow" interaction: VM(U Q , U p ; 0)+VN(U p , U Q ; 0). Although it does 
not have the same cancellation property as "fast-fast" interaction, the next lemma reveals that, 
upon averaging in time, it is of 0(e). In other words, "fast-slow" interaction is asymptotically 
negligible when averaged in time. 

Lemma 3.4. For any U = U Q + U p G X m with m > D/2 + 3 that solves $tt$ and satisfies 
the non-vacuum condition ||e/o||c([o,T];_f/ m ) — 1/2? 

rT 

VM(U P , VU Q ; 0) + VM(U Q , VU p ;0) dt 

Hm-2 



< Cmax\\U\\ H m yL + T max \\U\\h-" 
Proof. Since the definition of V is independent of the density component, it suffices to estimate 



V 





l p -Vu q + u Q -Vu p 



dt. 



By (|2.10p . V can be dropped as well. Now, we demonstrate the estimate on u^-Vu p , the other 
part being very similar. 

The idea is to treat the slow component u p like a constant and focus on averaging the 
fast oscillatory part in time. Integrating by parts is therefore in order. But first, a quick 
calculation reveals that 



by CV = 
by HSU 



CU Q ds= I CUds 
Jo 

e [ d t U + M(U,VU;e)ds 
Jo 



Since Proposition[23]yields \\U Q \\ Hm < C\\£U Q \\ Hm -i and therefore || /* U Q \\ H m. < C\\ jj £U Q \\ H m-i, 
it follows 



XjQ 



< Ce 



U 



+ / M(U,VU;e)ds 



H" 



By ([T5]> . the RHS is bounded by 

... < Ce (max ||£/|| m _i + imax ||[7||^m (||J7||^m + \\g(p; £)\\h™) 



\ t 
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The definition of g in (j3.3j) . together with p'(l) = 1, indicates that its dependence on e is 
essentially via positive powers of ep. The non-vacuum condition ||ep||c([o,T];_H" m ) — 1/2 then 
ensures that ||#(p; e^lg-m-i < C||p||^m-i. Therefore, it has been established that 

ft 



jjQ 



H" 



< Ce |^max ||f7||#m + t max ||C/||j^ n 



and, in particular, 



:= / v,® is bounded in H m (£l) norm by the RHS above. 
Jo 



Then, perform integrating by parts in time as follows, 

T T 

u Q Vu p dt = ii Q Vu p - 



and arrive at 



H r ' 



u Q Vu p dt 



u Q -Vd t u p dt, 



<Ce |max ||?7||H m + Tmax 



^max ||{7 p ||/fm + T max\\dtU p \\ Hm -i 



The last term of ||^C/ p ||^m-i are estimated using (|3.6p and the fact that M depends on e via 
positive powers of ep. □ 

To close this section, we prove Theorem 13.11 

Proof of Theorem \3.1[ Lemma 13.31 and 13.41 together with the integral version of (|3.6p lead to, 



for — e(t, •) :=u i 



+ 



T(u p -Vu p )ds, 



<Cemax ||C/||^m ( 1 + Tmax 



(3.7) 

'■" i \ 

Here and below, V is merely the restriction of V on the velocity component. In differential form, 
e satisfies 

(3.9) -d t e = d t u p + V(u p -Vu p ), e(0, •) = (). 

Then, subtract (|3.9p from the incompressible Euler equations (|3.4p and use the fact that I—V 
always yields a potential flow to obtain an equation for e := u — u p , 

d t e =d t e + (e + u p ) -V(e + u p ) - u p -Vu p + Vq, 

where q is determined from the incompressibility condition V-e = V-e = 0. 

Directly applying the energy method to estimate e in the above equation won't work since 
dte is of order 0(1) by (13. 9p . To this end, we "hide" dte and look at e\ := e — e which satisfies 

-d t ei{t) =(ei + e + n p )-V(ei + e + u p ) -u p -Vu p + Vq 

=ei -Vei + e x -V/ + /-Vei + f + Vq. 

Here, / := e + u p , f e := (e + u p )-V(e + u p ) — u p Vu p . By (13. 8p . they satisfy estimates 



(3.10) 



< i^(max T), 



\\f £ \\ H m-s < sF(max\\U\\ Hm ,T) 
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for some polynomial F. 

To apply the energy method to the above system, note that e, e and e\ all belong to H m ~ 2 (Q) 
and satisfy the incompressibility condition and solid-wall boundary condition. Thus, by elliptic 
estimate (|2.ip . 

lleiHfl-m-8 < C(||ei|| L 2 + ||Vxei|| ff m-4). 

It remains to establish upper bounds for dt ||ei \\i2 and <9t||Vxei || £J m — 4 111 terms of ||£7||iim- For 
simplicity, we skip this part, only alluding to the facts that (|3.10p is a symmetric Hyperbolic 
PDE system with dissipative boundary condition and is equipped with a vorticity equation 
for Vxei. Also, one can take spatial derivatives on (|3.10p up to (m — 3)-th order and still 
retain H l (£l) regularity for <9™ _3 ei, d™~ 3 f and therefore H 1 l 2 (dVL) regularity for any boundary 
integrals involved in the process. 

Finally, we obtain estimates for u — u p = e\ + e as desired. □ 

4. A Priori Energy Estimates 

Recall the equation of state: pressure=p(l + ep) with p(l) = 0, p'(l) = 1 and introduce a new 
unknown r := p ( 1+£p ) g.t. 



e 



1 / „ s V-u 

(r t + u-Vr) + = 0, 



p 1 (er)p' (p~ l (er)) 



(4.1) 

Then, rescale this system by replacing 



_i Vr 
p (er) (u t + u-Vu) H = 0. 



1 

u — > eu, r — >• er, t — > —t 

e 

and arrive at 

p^ 1 (r)p' (p~ l (r)) ^ 1 ^ ' 

p -1 (r) (u t + m-Vu) + Vr = 0. 

Finally, rewrite it as a symmetric hyperbolic PDE system for the rescaled variable V := [ 

\u 

D 

(4.2) A (V)d t V + Aj(V)d Xj V + £ M = 

j=i 

with the obvious definitions of Aq, Aj. In particular, ^4o(0) = ^d+i- 

This rescaled system will be the main subject of this section for it offers the convenience of 
free of any e terms. 

Define an energy norm for scalar or vectorial function f(t,x) at any given t, 



(4.3) 



I \d?dgf(t,x)\ 2 dx. 



Also, introduce some shorthand notations for simplicity. Let B$(t,x), Bi(t,x), ... Bj)(t, x) de- 
note symmetric-matrix-valued functions of dimension (D+l) x (D+l). Define B := (Bq, Bj}) 
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and diuf )X B := dtB^ + Y^j=i^XjBj. Define the boundary matrix using the outward normal 
n = (ni,...,n D ), 

D 

B n := n jBj on Sf2. 

3=1 

The main theorem of this section is regarding a priori estimates of linear systems equipped 
with certain structural features. Estimates for nonlinear systems then follow naturally. 

Theorem 4.1. Consider the linear symmetric hyperbolic system, 

D 

(4.4) B d t V + Bj^V + = ^ 

i=i 

(4.5) Han = °> 

with possible dissipation B n \ d ^ > and all the eigenvalues of Bo located in [1/2, 2]. Also, impose 
a structural assumption that there exists a "vorticity" operator K. of first order differentiation 
s.t. applying K to this system results in a "vorticity equation" (at least for smooth solutions) 

D 

(4.6) B Q d t K[V] + Bjd Xj K[V] = F + F x 

i=i 

where F\ is a sum of products of entries from c^B with entries from dt, x V. 

Then, for any V, B G Y!j=o C j ([0, T); H m - j (VL)) with integer m > D/2 + 1 satisfying (jOjh 
(|4.5p . there exists a universal constant C = C{m,Q) s.t. 

(4.7) \\V\\ m {t) < c{\V Q \\ m + \\F\\ m + S{B)\\V\\ m + J £{B)\\V\\ m + \\F\\ m dt') 



where 



D 

£{B) := \ld t>x B \l m -i + YM 



^ lll-^jlllm- 
J'=l 

Remark 4.2. The compressible Euler equations naturally satisfy the structural assumption due 
to the vorticity equation 

dfio + u-Vu + lower order terms = 0. 

This theorem is an immediate consequence of the following two lemmas. Lemma 14.41 uses 
energy methods and a careful mollification strategy to estimate the space-time norms of the 
vorticity K.V and the 1? norms of d\V . All boundary integrals in this calculation either vanish 
or have the correct signs since: 1. the vorticity equation does not have an C term; 2. time 
differentiation of V retains the boundary condition, at least up to (m — l)-th order. Lemma 
14.51 uses elliptic estimates to bound the rest part of |||V||| m in terms of the already established 
bounds. 

To accommodate the lack of boundary regularity for d+^V G L 2 (Q), a mollification process is 
used in the interior of We note by passing that it is proved in Rauch [18]: V being a weak 
L 2 solution to certain type of hyperbolic PDE system implies the trace V-B n V exists in a very 
weak sense. For the mere sake of a priori estimates, nevertheless, we provide a self-contained 
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argument with somewhat loosened assumptions. To this end, pick a C°° function ip(-) with 
suppip = (—1, 1), < ip(-) < 1 and J R ip{z) dz = 1. Then, define a C°° mollifier in the x variable 
ipg(x) := nJLi ^ (x) anc ^ an °ther one in the t variable ip^t) := ^ip (|). 
Also, define fi^ := {x € $7|disi(x, 9f2) > 5}. 

Proposition 4.3. Let x\{t,x) := 1[s,t-S]xU and XfCM) '■= ^O.TlxfV Then, for any f{t,x) € 
C([0,T];L 2 (O)), as 6 ^ + , 

xi-m*f) 

both converge to f(t,x) strongly in L 2 (Q,) for fixed t € (0,T) and strongly in L 2 ([0,T] x $7). 
Moreover, for any C 1 function g(t,x), the commutators 

xl-dt {g4 *f-4* (gf)) , xs-9 x (#f * / - r s * (gf)) 

converge to weakly in L 2 ([0,T] x Q). 

The proof is classic. We remark that the weak convergence is a direct consequence of the 
strong convergence since one can always choose a testing function in Cq([0,T] x Q) to allow 
integration by parts. This suffices to validate weak convergence since Cq is dense in the L 2 
topology. 

Lemma 4.4. Under the same assumptions as Theorem\4-l\ 



t2 



* 2 / r* 2 



(4.8) l\JCV\l (m _ 1:Bo) ^ + |||d f m Vl(o,B ) ^ < C (J^ \ld t , x B\\\ m ^l\V\l m + |||F||| m dtj , 

for < t\ < t2 < T. On the LHS, the subscript Bq is understood as a weight that modifies the 
vector norm | • | in the definition (|4.3p . 



Proof. By continuity in time, we only need to consider < t\ < t% < T. Also, it is easy to check 
that for any continuous and positive-valued function f(t), g(t), h(t), 

](t)f 2 (t) + h(t)f(t)dt. 



f\x < r g{t)f{t)+h{ t )dt iff f%\ <2 r 9 { 
j t± j ti 



Thus, it suffices to prove 

II 2 



(4.9) ||x;^ m _ 1)Bo) ■ ' ' 2 f /,J 



I + \Wvf %Bo) 



ti 



<C[ I Ilia ^IIL-! Illy III 2 + |||F||U||y||| m d* 



To estimate JCV, first apply df x with |a| < m — 1 to (|4.4p and obtain a system for V a := df x V 
well defined in C([0, T]; L 2 (S7)), ' 

D 

(4.10) B d t V a + Bjd X] V a + C[V a ] = F a 

i=i 

where 

d H 

3=0 \P\=1 
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with eg = if a — j3 contains negative index. Next, convolve this L 2 system with rpg in x variable 
and apply /C to obtain, 



D 



(4.11) 



B d t w + Bjd Xj w = Gx + G 2 



where 



w :=/C^| * V c 



G- 2 



D 



V° 



for (t,x) G [0, T] x Qg. Here, we used the facts that ICC = and everything is sufficiently 
smooth. 

Now, take L 2 ([ti,t2] x Qg) inner product of (|4.1ip with w, apply the divergence theorem in 
t, x and use the symmetry of Bq, Bj) to calculate, 



/ w-(Gi + G 2 )= / wB d t w + w-Bjd Xj w 

J[ti,t 2 ]xn s Jlh^xQg j =l 

I (w-B d t w + y^w-Bjd Xj w + w-(div t . x B)w) + [ w-B w + J 

JlhMx^s^ j= i J J^s h J[t 



'[*i,t2]xn a 
It follows that 



+ / w-B n w 

h 7[ti,t 2 ]xan« 



(4.12) 



W-BqW 



t-2 



2w-{G\ + G2) + w-(div ttX T5)w J — / w-B n w 



Here, outward normal n is extended smoothly from dSl to f2 s.t. it is also the outward normal 
of dVtg for 5 < 1 /curvature of <9f2 and vanishes faraway from d£l. 

Apply Proposition 14.31 on (|4.12|) to pass the strong limit as 5 — > + on both sides. More 
precisely, up to a multiplication of Xg > 

w(t,-) -»• d a (ICV{t,-)) strongly in L 2 {£1) for fixed t 



w -»■ <9 a (/CF), Gi -> /CF° strongly in L 2 ([0,T] x 0) 



and therefore the LHS of ()4.12j) and the integral of w-G% and w-divt, x Bw on the RHS converge 
to quantities that are bounded by corresponding terms of ()4.9p . 
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Now, we are left with the integral of w-G2 and the boundary integral. Subtract and add 
JC(Bodt + J2j=i Bjd Xj )ijjg * V a in between the two commuting terms of G2, 

D D 

G 2 =(B d t + Bjd Xj )W x s * V a - Kitf * (B d t + ]T Bjd Xj )V a 

3=1 3=1 

D D 

= (B d t + Bjd Xj )JCi> x * V a - JC(B d t + ]T Bjd Xj )i>$ * V a 

D D 

+ lc((B d t + Bfi^Ws * V a ~ V>? * (B d t + B 3 d X] )V a ) 
3=1 j=i 
= :h-I 2 +h-h- 

By the structural assumption (|4,6p . I1 — I2 consists of a sum of products of entries from c^B with 
entries from dt^ips * ana - therefore is bounded as desired. The I2 — I3 part, by Proposition 

31 converges weakly to in L 2 ([0, T] x Q). We then have shown lim sup / w-G2 dxdt is bounded 

5^0+ J 

by the RHS of flU}. 

And lastly, for the boundary integral — ju t2 ] x gn s w-B n w, since n is chosen to be the outward 
normal of dtls for < 5 < l/curvature(9r2), there exists a straight path Xq := xq + 5 n(xo) for 
< 8 < 1 that connects any point xq = Xq £ dQ$ to the closest neighbor x\ = X\ £ d£l. Some 
Geometric Analysis reveals that n(xo)|g 05 = n(xi)| 9f2 . It then follows 



/ —B^xie^de^cs^BU 

Jo <9 n 



B„(xi) - B n (x ) = (Xi - Xq) 

Then, together with the assumption B n \ dQ > 0, it guarantees 

(4.13) - [ w-B n w <C5 T IQeBloolMl^an 

J[ti,t 2 ]xdn s Jt! 



'[ti,t 2 ]xdn s 

Furthermore, the Divergence Theorem implies 



H\h(n s ) 



/ (|w| 2 n)-n=/ V-(\w\ 2 n) <C\\w\\ L 2 {ns) \\w\\ H i {ns) 



Use the definition of w and the scaling invariance IjV^f = C/5 to continue, 
... < C\\KV a \\ L 2{l + ||VV?|| Ll )l!/Cm L2 < (C/5)\\V\f m . 

Combined with ()4.13j) . it shows the boundary integral is bounded as desired. We have finished 
showing the estimates on ICV in (|4.9p . 

The study of S^~ X V involves a very similar strategy. Let a = (m — 1,0, ....0) in (|4.10p and 
convolve it with -^ipg, 

D 

(4.14) B d t w + Bjd Xj w + £[w] = G l + G 2 

3=1 
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where 



w --dtiPl * d™~ l V = d^ips * V 
G x :=d t ^ 5 * F a 



G 



2 ■ = 



D 



3=1 

This system is well defined for (t, x) £ [8, T — 8] x £1 and, in particular, w retains the same 
regularity as V. Thus, w is C 1 and satisfies the same solid-wall boundary condition 

0. 



w-n\ 



an 



Now, follow the same procedure that leads to (|4.12p . knowing that it can be performed all the 
way up to [8, T — 8] x £2 and the C\w\ term in (|4.14j) disappears as J n w-C[w]dx = due the 
above boundary condition, 

' / 2w\Gi + G 2 ) + w\div t ^ x B)w 



(4.15) 



W-BqW 



*o J[to,ti]xQ 
< 



w-B n w 



[to,ti]xan 



2w-(G x + G 2 ) + w-(div t , x B)w _ 

' [to,ti] xfi 
for 5 < t < ti < T - 5. 

Fix t\,t2 and pass the limit as 8 — > + in (|4.15p just like what we did for (|4.12j) . we prove 
AMD for < h < t 2 < T. 

□ 



Lemma 4.5. Under the same assumption as Theorem \4-l\ 

\\d^v\\ m _ k < c(ix:v|| m _i + ||a f m ^|||o+ll^|m-i + 11^1 

(4.16) 



D 



+ (II^Bo| m _i+5^|lS J -||| m )|||F|| 

3=1 

for any < k < m. 

Proof. We perform induction w.r.t. k. 
First, when k = m, it is trivial. 
Assume (|4.16p is true for k = n + 1 

(4.17) iar +1 y| m -n-i < RHSdns}. 

Apply <9 t n on (IQjl . 

(4.18) B d? +1 V + C[d?V] = F' 
where 

D 

F' := F — d?{B Q d t V) + B d? +1 V + d? B 3 d X] V) . 

3=1 

It is easy to check that |||.F / ||| m _ n _i < RHS (|4.16|) by virtue of (|1.8|) . Likewise 
iBo^VlU-n-i < iBolool^rVlU-n-i+RHSdllS]) 
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Due to the inductive assumption (|4.17j) and the assumption 1/2 < Bq < 2, this term is also 
bounded as desired. Plugging all these estimates into f|4. 18[) yields 

(4.19) \\C[d?V]i m - n „i < RHSgH]). 

Finally, note that elliptic estimate (|2.4p implies, 

|0*(a?tOIIm~n-i <c(i\c[d?v] | ro _„_ a + \iJC[d^v]\i m ^x + |^yi| m _ n _i). 

Combined with (|4.19p . it yields the bound on ||3a;(9"V)|| m _ n _i as we expect. Since all the 
highest order derivatives in |||9"V||| m _ n contain at least one x derivative unless it's d™V, the 
above estimate sufficiently leads to (|4.16p for k = n. □ 

Obviously, Lemma 14.41 and 14.51 also apply to any lower order derivatives that appear in the 
definition of |||V||| m . The very lowest norm |||V|||o is merely L 2 and can be bounded using classical 
approaches. Therefore, combine all these estimates and finish proving Theorem 14.11 

Now, for the nonlinear system (fO) . we set B = (A (V), Ai(V), Ad(V)) and F = in 
Theorem 14.11 All conditions are satisfied, esp. the structural assumption (|4.6p due to the 
existence of an actual vorticity equation for the compressible Euler equations. We then use the 
definitions of £ and Aj(V) = Ao(V)uj to arrive at: there exists a polynomial 7r(-) of degree 
m — 1 with positive constant coefficients only depending on m s.t. 



^(i?)<^l4 fc) (^)lool^lUvr(|y| 



k=0 

Since Aq(0) = / and is defined using the pressure law p(-) G C°° with p(l) = 0, p'(l) = 1, there 
exist constants C\,Ci dependent on p(-) and m s.t. 

m 

^l4 fc) 00ioo < Ci and l/2<A (y)<2 if [^U < C 2 . 

fc=0 

Plug it into (|4.7p to arrive at, for any |||V||| m < C 2 

tv\\ m (t) < c(|||Fo|||m + Ci7r(a 2 )|||y|||^ + J*cMc 2 wf m di/) 

which further implies, 

1 /** 

\lV\l m (t)<2C\\V \l m + — / illT^IH^ cZ*' 

^3 Jo 

if |V| m is also bounded by C 3 = 1/{2CC 1 tt(C 2 )), 

By a simple comparison argument, one can show that |||V||| m (f) < a(t) for a(t) solving 

1 f* 

a(t) = o(0) + — / a 2 (t') dt', a(0) = (2C + l)j||Vb||U 

^3 70 

as long as |||V||| m (i) < min{C2, C3}. Since this integral equation of a(t) is equivalent to a Riccati 
type ODE, we solve to get a(t) = t 1 t . Thus, 

|||V|||m(i) < j — as long as the RHS < min{C2, C3} 

(2C+l)|||Vb||| m ~ 



LOW-MACH-NUMBER EULER EQUATIONS 



19 



Finally, since all the time derivatives in the definition of |^|m can be expressed in terms of 
spatial derivatives via system (14.21) . one has 



II II m < C4 1 1 "Vol | if™ if ||Vo||// m < Cs- 

Theorem 4.6. Let m > D/2 + 1. There exists constant C* , C** , C*** only dependent on 
m,Q and the pressure law s.t the rescaled compressible Euler equations (|4.2|) with initial data Vq 
satisfying the compatibility condition admits a unique solution in the class of 



E c '( °<Ff^— ■,H m - j (n)) if \\v \\ Hm <c**. 

n VL \\vo \\h™ ' 
3=0 



The solution is uniformly bounded \\V\\ H m^ < C***\\Vo\\H m on this finite time interval. 

Consequently, the fast oscillatory system (14. lj) with initial data (ro,Uo) satisfying the compat- 
ibility condition admits a unique solution in the class of 



C* 



E ci (k in— l Hm - j W) if \\(ro,uo)\\ Hm <C**/e. 

The solution is uniformly bounded \\(r, u)\\H m (Q) — C***ll( r 0) Wo)||_f/ m on this finite time interval. 

Here, we assumed the existence theory is automatically valid for the compressible Euler sys- 
tem. For linear Hyperbolic PDE systems with characteristic boundary, the existence is studied 
in Lax & Phillips [15] in terms of strong L 2 solutions and in Rauch [18] for H m regularity. The 
nonlinear case is studied in Schochet |20U21| . Our Theorem 14 . 1 1 14TB1 and their proofs provide clear 
evidence that the continuation principle can be performed on a uniform time interval indepen- 
dent of e. In particular, when a <5n-Vu term is added to convert the original singular boundary 
matrix to a nonsingular one, the uniform estimate in Theorem 14.11 remains valid because: 1. 
boundary matrix is still dissipative; 2. a vorticity equation still exists. 

We also note by passing that a crucial "maximal positivity" condition on the boundary matrix 
and boundary condition is needed for all these existence theories to work (cf. Rauch [18] and 
Schochet [21]), but the system dealt with here is a canonical case satisfying this condition for 
sufficiently small e. 

5. Extension to the Rotating Shallow Water Equations 

We can extend the above framework to the 2D Rotating Shallow Water (RSW) equations 
in a very natural way. Consider a thin layer of fluid moving under the gravitational force and 
Coriolis force. Assume the fluid is vertically homogeneous so the dynamics is essentially in 

2D. Let u = ( Ul ) be the velocity field and p be the height perturbation against a constant 

\U2j 

background height. Then, we formulate 

(5.1) d t p + V-(fm) + -V-u = 0, 

(5.2) d t u + u-Vu + ^Vp+^u ± = 0. 
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Here, u 1 - := [ u<2 \ reflects the Coriolis effect due to a rotating frame. Impose boundary 
V u x ) 

conditions as 

u-n| =0, / p = Q. 

Jn 

The analogue between the above system and the compressible Euler equations is as follows. 
The pressure law for the RSW equations is p(p) = \p 2 , which is due to the gravitational force. 
The singular parameter e is the 2D Froude number and plays the same role as the Mach number. 
The elliptic operators, with the same notations as before, are 



CU 



Vp + u A 



KJJ := d x ii2 — dyUi — p. 
The RSW system is endowed with a vorticity equation 

(5.3) d t (ICU) + V-(plCU) =0. 

It remains to justify the Propositions in Section [2] regarding elliptic estimates and orthogonal 
projections. The elliptic estimates part is quite straightforward since the C, IC defined here 
are only lower order perturbation of their counterparts in Section [2j and one can always bound 
norms of p by norms of Vp. For the orthogonal projection part, we essentially need to find 
an elliptic PDE to define V, Q which is as convenient as ()2.7p . Indeed, since (p,u) £ Ker£ iff 
u = (Vp)^, we characterize V in the following way, 



VU = I / v ^± j with ( A - = d xU2 ~ d. y ui - p, 4>\ dn = 0, 
and, in a much better form, 

(5.4) V = K.*(JOC*)- 1 /C 

where IC* is the formal adjoint of /C and subscript D indicates the Dirichlet boundary condition 
associated with the elliptic inverse. 

The identity (|5.4[) suggests a parallel argument in Linear Algebra and, in particular, the 
method of Least Squares. Without proofs, we state the following set of duality relations which, 
unlike the Euler equations, are independent of domain topology, 

ImgQ = Ker/C = Img/3* =KerP, 

IrngT 7 = Ker£ = Img/C* =KerQ. 

Here, C* = -£. 

The last point we would like to make is the analogue to Lemma 13.31 regarding cancellation 
of "fast-fast" interaction under the projection V. This is still valid for RSW equations due to 
a simple observation: the vorticity equation (|5.3p results from applying IC to (|5.ip . (|5.2p and 
using ICC = 0. In other words, the identity 
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is true regardless of the equation. Now, let U = be the fast part of a solution. Then, by the 
duality relations, fCU® = 0, and the above identity becomes 



Again, by the duality relations, we obtain 
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